In-Plane Dynamics of Membranes with Immobile Inclusions 
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Cell membranes are anchored to the cytoskeleton via immobile inclusions. We investigate the 
effect of such anchors on the in-plane dynamics of a fluid membrane and mobile inclusions (proteins) 
embedded in it. The immobile particles lead to a decreased diffusion coefficient of mobile ones and 
suppress the correlated diffusion of particle pairs. Due to the long-range, quasi-two-dimensional 
nature of membrane flows, these effects become significant at a low area fraction (below one percent) 
of immobile inclusions. 

PACS numbers: 87.16. dj 87.14.ep 47.15.gm 47.56.+r 



Biomembranes are key structural and functional in- 
gredients of any living cell. They are based on a self- 
assembled fluid bilayer of amphiphilic molecules (lipids) , 
containing a high concentration of embedded proteins 
that perform vital biological tasks [l|. From a physical 
viewpoint, biomembranes represent an interesting state 
of matter, whose properties are intermediate between 
two-dimensional (2D) and three-dimensional (3D), and 
between fluid and solid. In the past four decades the 
structural properties of membranes and the dynamics of 
their out-of-plane fluctuations have been well character- 
ized [2|, |3|. The in-plane fluid dynamics, which is cru- 
cial, e.g., for the diffusion of membrane proteins, was 
first addressed by Saffman and Delbruck (SD) Their 
model considered the motion of a single inclusion within 
a viscous slab surrounded by an infinite fluid. While cer- 
tain experiments were consistent with the SD predictions 
other works pointed to discrepancies The SD 

model was extended in various directions, e.g., to treat 
proteins of larger size [9[ and number [lfl], viscoelastic 



effects [ill, [l2j, in- plane concentration fluctuations and 
phase separation I13l4l6l ] . and the effect of an adjacent 
rigid surface 17H13|. Contrary to these models, most 



biomembranes are neither freely suspended in an infinite 
fluid nor attached to a rigid wall. Rather, they are an- 
chored via proteins or protein-associated domains to a 
dilute, soft network of filaments (cortical actin in ani- 
mal cells or cortical microtubules in plant cells) [l[ . Our 
aim is to explore how the membranal in-plane response, 
and the consequent dynamics of mobile inclusions, are 
affected by such immobile anchors. 

Models of membrane dynamics are divided into those 
that conserve momentum in 3D 

USB! EMI, as appli- 
cable to freely suspended membranes, and those that do 
not 17 - as appropriate for substrate-supported mem- 
branes. This essential difference is reflected in the respec- 
tive velocity Green's functions, G froc (r) and G sup (r)- 
tensors that give the flow velocity of the membranal fluid 
at position r, v(r), in response to a localized in-plane 
force, F, applied at the origin, Vi (r) = G l :; c/sup (r)F 3 
(with i,j = x, y and summation over repeated indices). 
In Fourier space, G(q) = J d 2 re _lqr G(r), these are 
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where we have assumed a flat, incompressible membrane 
with 2D viscosity ?7 m . In Eq. ([1]) k _1 is the SD length, 
the characteristic distance beyond which the free mem- 
brane exchanges momentum with the surrounding fluid 

In the SD model kT 1 = rj m /(2r)f), where rji is the 
viscosity of the outer fluid. Typical values for the mem- 
brane 2D viscosity are rjm ~ 10 _10 -10 -9 Pa-s-m, yielding 

~ 0.1-1 /im. In Eq. (J5J) oT 1 is the distance beyond 
which the supported membrane loses momentum to the 
substrate. (For a membrane lying a distance h from a 
rigid wall, a = [k/(2/i)] 1 / 2 [Hj].) Despite the appearance 
of the decay parameters n and a, Eqs. {J) and ([2]) de- 
scribe long-ranged velocity responses — G free decays as 
1/r and G sup as 1 jr 2 (due to the conservation of 3D mo- 
mentum in the former, and 2D membrane mass in the 
latter [Hi). 

Biomembranes contain anchoring inclusions separated 
by typical distances of 30-80 nm [H, i.e., covering an 
area fraction <j> ~ 10 _3 -10 -2 . These immobile inclu- 
sions break the membrane's translational symmetry and 
absorb momentum; hence, the large-distance response 
should be similar to G sup . From another perspective, 
we expect the immobile obstacles to serve as a porous 
matrix within which the lipids flow. Indeed, Eq. ^ is 
analogous to the response of a fluid embedded in a porous 
medium 23j. We establish below that this intuitive pic- 
ture is correct and derive the momentum decay length 
as a function of the area fraction qb of immobile inclu- 
sions. As this length decreases from arbitrarily large to 
smaller values upon increasing <p, the membrane crosses 
over from a dominantly "free" behavior to a "supported" 
(or "porous") one. 

We first note the quasi-2D nature of the system. Over 
sufficiently small distances a membrane behaves as a 
momentum-conserving 2D fluid, whose velocity Green's 
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function is given by 24 1 
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where /3 _1 is a cutoff length regularizing the divergent 2D 
behavior and 7 ~ 0.577 is Euler's constant. Both Eqs. ([1]) 
and @ , upon inversion to real space and taking the limit 
of small distances, coincide with G 2D , with j3 = k or j3 — 
a, correspondingly. In the current problem the immobile 
inclusions produce an effective cutoff for the membrane's 
2D behavior, to be derived below. With increasing <fi, 
P~ l decreases from k~ 1 down to a microscopic length. 

Let us consider a bare free membrane and ask how its 
velocity response, Eq. ([T]), is modified by the presence 
of immobile cylindrical inclusions [25| . The inclusion ra- 
dius, a, is taken as the smallest length in the problem, 
and the calculation is restricted to the leading terms in 
Ka and a/r. We begin by applying a localized force F 
at the origin. The resulting flow velocity of the bare 
membrane is v!f\r) = G^ e (Y)Fj. Next we place an im- 
mobile inclusion at position r'. Due to the flow v^ ', the 
inclusion exerts an additional force on the membrane, 
F« ~ -rvW(r'), where 



r = 47rW{ln[2/(/?a)] - 7} 



(4) 



is the friction coefficient of the cylindrical [25[ inclu- 
sion @, [H HI, /3 being an effective 2D cutoff. We 
neglect moments of the force distribution higher than 
this monopole. The inclusion-induced force changes the 



flow velocity by vf\r) = Gff c (r - r')F^'(r'). Thus, 
a single immobile inclusion reflects the original force as 



.(1), 



i^(r) = -TGY° c (r~r')GJ£ c (r')F k . Now, consider many 
randomly distributed immobile inclusions, covering an 
area fraction d> of the membrane. To first order in <b, their 



effect on the membrane's flow is (v\ (r)) = G\j(r)Fj, 
with G«(r) = -r[0/(7ra 2 )]/d 2 r'Gfr(r-r')Gl7 c (r'). 



At a higher value of (f> many-body terms set in — the 
flow reflected from one inclusion is reflected again from 
another, and so on. We have calculated all orders of these 
monopolar hydrodynamic terms while continuing to as- 
sume a uniform static distribution of inclusions [26j . This 
yields fj(r) = Gfj(r)Fj, with a cumbersome expression 
for G eff (r) HI. In Fourier space, 



(q) = 
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Vm[q(q + k) + A 2 ] 

A = a [r0/(^ m )]- 1 / 2 . 
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In the limit r 0, G has the expected form of G 
[Eq. (|3j)]. with a cutoff /3 that satisfies 
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Equations (jlj), ©, and ([7} provide a self-consistent 
scheme for obtaining /3a and A/a as functions of <f> and 



Ka. The decrease of /3 _1 and A with increasing is shown 
in Fig. [1] In the limit cf> we have A — > 00 [Eq. (J6j)] , 
and Eq. ([7]) then gives f3 —> k (see Fig. [TJ\) — i.e., the 
cutoff of the 2D behavior is that of a free membrane, as 
expected. As 4> increases, A decreases roughly as -1 / 2 
(Fig. [TJ3 inset), while /3 _1 decreases initially as (/)ln</>, 
and subsequently more sharply (Fig. [TJA.) . Both lengths 
reach values comparable to the inclusion size a, where- 
upon the theory breaks down. The breakdown is marked 
by a loss of solutions for the self-consistent scheme, oc- 
curring at a fixed (up to corrections of order Ka), small 
area fraction, max ~ 0.058. This analysis reveals two 
distinct regimes (cf. Fig.[T]): a low-concentration regime, 



<§; (Ka) 2 1 ln(«;a) | , in which A ^ 



and a 



higher-concentration regime, (ko) 2 \ ln(«;a)| <C (f> < 
where A ~ f3~ 
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FIG. 1. Two-dimensional cutoff length (A) and momentum 
decay length (B) as a function of area fraction of immobile 
inclusions. Solid, dashed, and dotted curves correspond to 
Ka = 10~ 3 , 1CF 2 , and 10 _1 , respectively. The lengths are 
scaled by the inclusion size a. The inset of panel B shows A, 
scaled by the SD length on a log-log plot, demonstrating 
a roughly <^~ 1//2 decay. 

We can readily examine the effect of <\> on the 
self-diffusion coefficient of a mobile inclusion of ra- 
dius b. From Eq. and Einstein's relation, D s = 



3 



k B T{\n[2/(/3b)} - 7}/(47r?? ro ), where k B T is the thermal 
energy. Because of the increase of /3 with </>, D s decreases 
from its free value D f J cc as given by SD [ij (the expression 
above with f3 = n, yielding typical values of Dl TQC ~ 1-10 
/um 2 /s), to significantly lower values (Fig. [5]). 
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FIG. 2. Self-diffusion coefficient of a mobile inclusion as a 
function of area fraction of immobile inclusions. The coeffi- 
cient is scaled by its value in a free membrane (0 = 0). The 
radii of the two types of inclusion are taken as equal, b — a. 
Solid, dashed, and dotted curves correspond to na — 10~ 3 , 
10 -2 , and 10 _1 , respectively. 

Next we address the coupling diffusion coefficients 
of a pair of mobile inclusions separated by a distance 
r. The longitudinal diffusion coefficient, £>l(0 = 
(Aa;iAx2)/(2f), and the transverse one, D^{r) = 
(AyiAj/ 2 )/(2t), characterize the correlated diffusion of 
the pair, respectively, along and perpendicular to their 
connecting line. Here Axi and Ay^ are the displacements 
of particle i along and perpendicular to the connecting 
line during time t. In the limit r» t, the pair diffusion 
coefficients are directly obtained as -DlM — ^BTG^(rx) 
and -Dt(t') — fceTG^rx), where x is a unit vector along 
the connecting line. As is evident from Eq. ((5|), the cou- 
pling diffusion coefficients in this limit depend only on A 
and k; they are independent of the size of mobile inclu- 
sions and depend on the size of the immobile ones only 
indirectly, through A. 

The results obtained for D^(r) and Df{r) using the 
full expression for G eff (r) (2(| are shown in Fig. [3] (solid 
lines). They include several simpler asymptotic regions. 
The low-concentration regime (kA 3> 1, Fig. |3l\) in- 
cludes three such regions. At short distances (r <C ft" 1 ) 
the behavior is 2D-like, resulting in D^,Dt ~ |ln(/3r)| 
with /3 ~ k (dashed lines). At intermediate distances 
(k _1 <C r -C A) the coupling becomes 3D-like, with 
Dl ~ 1/r and Z?t ~ l/?" 2 (dash-dotted lines). In 
these two regions the coefficients coincide with those in 
a free membrane, as obtained from Eq. (J]) [l(| ■ At suf- 
ficiently large interparticle distances (r 3> A) the cou- 
pling becomes sensitive to the immobile inclusions, with 



£>l, Df ~ ±l/r 2 (dotted curves in the inset). These co- 
efficients coincide with those in a supported membrane, 
as obtained from Eq. ([2} with a = A -1 [3]. The higher- 
concentration regime (kA <C 1, Fig. |3j3) contains two 
asymptotic regions. At short distances (r -C A) the be- 
havior is again 2D-like, Z?l,-Dt ~ |l n (/?''')|j but with 
f3 ~ A -1 S> k (dashed lines). At large distances (r S> A) 
we have again D^ 7 Dt ~ ±l/r 2 (dotted curves) as a re- 
sult of momentum loss to the immobile inclusions. 
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FIG. 3. Coupling diffusion coefficients of a pair of mobile 
inclusions as a function of mutual distance. Panels A and B 
present, respectively, results for low concentration (kA = 100) 
and high concentration (kA = 0.01) of immobile inclusions. 
Results obtained from the full Green's function (solid lines) 
are shown together with their various asymptotes (see text). 
The coupling diffusion coefficients are scaled by teT/i)m, and 
the distance by either the Saffman-Delbriick length k _1 or 
the momentum decay length A. 

Our predictions concerning the effect of immobile in- 
clusions on the self- and coupling diffusion coefficients of 
membrane proteins may be directly checked in particle- 
tracking experiments. The effective velocity Green's 
function derived above can be utilized [l4j in other the- 
ories involving anchored membranes. We have found a 
strong effect of immobile inclusions already at a very low 
area fraction (below one percent). The effect is manifest 
in the hydrodynamic screening length (Fig.[Tj3), the self- 
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diffusion coefficient of a mobile inclusion (Fig. [2]), and the 
dynamic couplings between two mobile inclusions (Fig. [3] 
A vs. B). The high sensitivity to the presence of immobile 
inclusions arises from the quasi-2D nature of the mem- 
brane and the resulting long-ranged flows. To highlight 
this special property, let us compare the situation to its 
3D counterpart. Consider a 3D fluid containing a volume 
fraction <f> of immobile particles of radius a. These ob- 
stacles introduce a momentum-decay length, A ~ a<f>^ x l 2 , 
which will strongly affect the self-diffusion coefficient of a 
mobile particle only when A becomes comparable to the 
particle size — i.e., for an appreciable value of <f>. Unlike 
the 3D case, which is characterized by a single (small) 
length scale, a, the quasi-2D dynamics depends on a 
and another, much larger length scale, /3 _1 . The long- 
range (logarithmic) nature of 2D flows [Eq. ©] makes 
this length sensitive to the immobile inclusions already 
for <f> ~ 1CP 3 (Fig. UK). The quasi-2D anomaly is re- 
flected also in the nonanalytic decrease of /3 _1 at small 
area fraction, as <f> l n <t>- 

The separation of length scales, (3~ x ^ a, has also al- 
lowed us to consistently include in the effective response 
of the membrane all the monopolar hydrodynamic terms 
(i.e., all orders of (f>) while continuing to ignore spatial 
correlations in the distribution of inclusions; for the low 
area fraction considered here (cj) < 10 -2 ) such spatial cor- 
relations are negligible. The other assumption necessary 
to make this expansion valid is that the distribution of 
inclusions is stationary, unaffected by the exerted forces. 
In addition, our analysis is restricted to the leading terms 
in na, a/r, and /3a. Corrections due to higher-order terms 
in the first two factors are related to higher moments of 
the force distribution exerted by individual inclusions. In 
actual membranes kcl is of order 10~ 3 , and the approx- 
imation is well justified. Higher-order terms in a/r be- 
come arbitrarily small at sufficiently large distances [13] ■ 
By contrast, since /3 is density-dependent, the restriction 
(3a <C 1 limits the validity of the analysis to cj) < 10~ 2 
almost irrespective of a (see Fig. UK). 

Immobile inclusions also affect the membrane's out-of- 
plane dynamics [28|, [2^] . To leading order in the normal 
deformations, the in-plane and out-of-plane dynamics are 
decoupled [111 ], and our quasi-2D model should remain 
valid within the deformed surface. At higher orders of 
the deformation, or because of inclusion-membrane cur- 
vature coupling, various corrections are expected 

We should mention two additional ingredients of ac- 
tual biomembranes. The first is the filamentous net- 
work to which the immobile inclusions are attached. 
Since the network is deformable, the anchoring inclusions 
should not be completely immobile. Qualitatively, this 
restricted freedom to move laterally will affect the mem- 
brane's in-plane dynamics as if the immobile inclusions 
had a somewhat larger effective size. In addition, the 
network adds a viscoelastic response, which is important 
for frequency-dependent properties but vanishes for the 



steady-state ones considered here. It also introduces an- 
other length scale — the network correlation length £ — 
which is of order 0.1 /im 1]. This large length scale will 
not significantly affect our results concerning the self- 
diffusion coefficient, but will influence the coupling coef- 
ficients at r > I. At such interparticle distances the im- 
mobile network can be replaced by an effective substrate 
[lij ]. The other ingredient is a finite density of mobile 
inclusions in addition to the immobile ones. Their effect 
can be readily incorporated by a proper renormalization 
of ?y m and k [lOj] . 

We have seen that for a reasonable value of <f> > 
(ko) 2 \ ln(Ket)| ~ 10~ 5 , the hydrodynamic interaction be- 
tween two mobile inclusions crosses over from a strong 
coupling [decaying only as ln(l/r)] for r < /3 , to a 
suppressed coupling (decaying as 1/r 2 ) for r > f3~ x (Fig. 
[3)3). We end with the following question. If the immobile 
inclusions were to divide the membrane into effective do- 
mains, such that within a domain two mobile inclusions 
are strongly coupled, whereas across domains the cou- 
pling is suppressed, what would be the value of <jp. Solv- 
ing the equation acf)^ 1 ^ 2 = /3 _1 5 we find <fi ~ 2 x 1CP 3 and 
5 x 10~ 5 for na = 1CP 2 and 1CP 3 , respectively. For a = 5 
nm this corresponds to domain sizes of about 0.1 and 0.7 
/im. These large values stem again from the quasi-2D 
nature of fluid membranes. 
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